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ABSTRACT 


Infrared geosynchronous satellite data with an interval of five 
minutes between images are used to estimate thunderstorm 
top ascent rates on two case study days. A mean vertical velocity 
of 3.5 ms"^ for 19 clouds is calculated at a height of 8.7km. This 
upward motion is representative of an area of approximately 10km 
on a side. Thunderstorm mass flux of approximately 2 x 10^^ gs"^ 
is calculated, which compares favorably with previous estimates. 
There is a significant difference in the mean calculated vertical ve- 
locity between elements associated with severe weather reports 
(w= 4.6ms'^) and those with no such reports (2.5ms'^). 

Calculations were made using a velocity profile for an axially 
symmetric jet to estimate the peak updraft velocity. For the larg- 
est observed w value of 7.8ms"^ the calculation indicates a peak 
updraft of approximately 50ms"T 

1. INTRODUCTION 


In the current paper we present the results of using a simple method to estimate 
thunderstorm cloud-top vertical velocity from SMS/GOES rapid-scan (5 minute inter- 
val) window channel infrared (IR) data. Time rate of change of cloud-top minimum 
equivalent blackbody temperature (Tgg) is converted to vertical velocity w by 


w = 



dt 


( 1 ) 


where the lapse rate is determined from rawinsonde data. In the following sections 
the calculated vertical velocities are compared for clouds with associated severe 
weather reports and for those with no such reports, and the computed values are also 
compared with previous estimates of thunderstorm vertical velocities. 

2. DESCRIPTION OF ANALYSIS TECHNIQUE 

The analysis of the digital satellite data was performed on the Atmospheric 
and Oceanic Information Processing System (AOIPS), an interactive image analysis 
system described by Billingsley (1976). Sequences of images are enhanced, ele- 
ments are isolated and identified, and their maximum gray level (minimum Tgg) 
recorded. The analysis of the thunderstorms used in this study (on April 24, 1975 
and May 6, 1975) was part of a larger effort (see Adler and Fenn, 1978). On each 
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of these study days an area of convection was monitored during a set period of 
time (approximately 4 hours). Thunderstorm elements were defined by the tech- 
nique described by Adler and Fenn (1976, 1978) and time histories of each ele- 
ment were determined. 

It should be emphasized here that not all thunderstorms can be observed at 
middle tropospheric heights. This is because they are often hidden by dense cirrus 
clouds produced by previous convection. In the two case studies that will be de- 
scribed here, of the thunderstorm elements defined above 10km (Tgg = 226 K), only 
about 25-30% could be detected at lower heights. 

3. VERTICAL VELOCITY ESTIMATES 

a. Sources of Error 

The vertical velocity estimates presented in this paper are subject to error be- 
cause of errors in the satellite radiance measurements and possible unrepresentative- 
ness of the data. A possible source of error is in the cloud emissivity. For thick 
water clouds such as we are dealing with, the emissivity (e) is very close to unity. 

In order to use Tgg in place of T in Equation (1), e = 1 must be assumed. For 
ice clouds, especially thin cirrus clouds, emissivity can be substantially less than 
1.0. As the thunderstorm top glaciates, the change from water to ice may reduce 
the cloud top emissivity and Tgg will be larger than T, the cloud top temperature. 
This effect, however, appears to be small and will not affect the vertical velocity 
calculations. This conclusion is based on calculations made with the help of tables 
presented by Hunt (1973) and on results presented by Cox (1977). 

The most serious source of error stems from the satellite data itself. The IR 
channel has an instantaneous field of view (IFOV) of 8km on a side at the sub- 
satellite point and approximately 10km at 40N. As noted by Negri, et al., (1976) 
the use of SMS/GOES IR temperatures to determine thunderstorm height results 
in underestimates, especially for small elements. In comparison with radar meas- 
urement of thunderstorm tops, the satellite underestimation is approximately 2 km. 
This effect is related to two factors. First, the satellite, because of its rather large 
IFOV, is averaging over an area approximately 100km ^ compared with the radar 
observation, which is applicable to an area closer to Ikm^. Thus the radar will be 
identifying smaller, higher features, because of its better resolution. The second 
factor is inadequate sensor response when going from a warm (low) to a cold (high) 
target (Negri, et al., 1976). In general, the bias in the estimation of storm heiglit 
will not affect the vertical velocity calculation. However, it may affect the height 
to which the velocities are assigned. 

Other errors might arise from the use of inaccurate lapse rates in Equation (1). 
In the calculations to follow we use a smooth profile which is a mean of the am- 
bient and the moist adiabatic lapse rate. The calculations can be shown to be 
rather insensitive to variations in the lapse rate. Using either the moist adiabatic 
or the ambient lapse rate instead of the average of the two produces only 10% dif- 
ferences in the calculated velocities. The final validation of the calculated vertical 
velocities must come through a careful comparison with radar, aircraft or satellite 
stereo observations. 

b. Vertical Velocity Results 

Fourteen elements on May 6, 1975 and ten on April 24, 1975 were analyzed. 
Minimum Tgg as a function of time for each element was plotted and dTgg/dt 
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values were calculated. The warmest or lowest of the monitored elements was at 
Tbb “ 260 K (about 6km). Not all thunderstorms could be observed from that 
point upward through the remainder of the troposphere. Some elements were 
blocked out by other storms; other elements were not detected until they pene- 
trated middle level cloud fields. 

Vertical velocities were calculated every 5K in the vertical for each cloud ele- 
ment or thunderstorm using Equation (1) and a lapse rate halfway between ambi- 
ent and moist adiabatic. The mean w was then calculated for various categories of 
clouds to produce composite profiles. The results for the May 6 case are shown 
in Figure 1 . Profiles are shown for thunderstorms associated with severe weather 
reports (based on National Severe Storms Forecast Center logs) and those with no 
accompanying reports. The numbers in parentheses indicate the number of cases 
constituting each mean or composite vertical velocity. The mean profile of all 
cases is also shown. 

The vertical velocities in Figure 1 and the calculated divergence noted above 
are applicable to an area of about 10km on a side. The calculated vertical veloci- 
ties do not represent updraft core velocities, which could be an order of magnitude 
larger when measured on a horizontal scale of I km (see Section 5). 

A similar diagram for April 24, 1975 case is given in Figure 2. The convection 
of interest on this day was centered in southwestern Missouri. The composite w 
profiles for three categories are displayed. The additional category is for weak ele- 
ments which did not reach a height of 10km (as determined by the Tgg values). 

In the layer from 7 to 9 km the average w is approximately 1.5 ms' ^ for those 
storms. This is significantly lower than the composite for “non-severe” elements 
in either Figure I or 2. 

The 235-240K level (^8.7 km) is representative of the layer of relatively large 
vertical velocities on both days. Figure 3 shows the frequency distribution of the 
19 elements or clouds for these two days. The hatched portion of the histogram 
contains the values for the severe weather elements. The average w for all 19 cases 
is 3.5 ms'^ The severe and non-severe elements have average values of 4.6 and 
2.5 ms'^ respectively. The severe thunderstorms dominate the high end of the dis- 
tribution where six out of seven cases with w > 4ms'^ are associated with severe 
weather reports. 

c. Mass Flux Calculation 

Because of the values of w calculated in the last section are representative of 
an area larger than a typical thunderstorm updraft, vertical volume or mass flux 
calculations can be simply made through the formula, 

= pAw, (2) 

where Fj^ is the vertical mass flux, p is the density, and A is the area. For the 
235-240K layer (~8.7km) p is assumed to be 5 x 10“^gcm'^, and A is assigned 
a value of lOOkm^ for the area of the satellite IFOV. 

With the given values for p and A the mean w for all storms of 3.5 ms'^ is 
converted to a mass flux of 1.8 x lO^^gs'^. The mean w of severe elements 
(4.6ms'^) is equivalent to a mass flux of 2.3 x 10 ^^gs"^. These magnitudes are 
for the mass flux through a given layer associated with a growing thunderstorm 
top. The calculated values compare favorably with results presented by other in- 
vestigators. Kropfli and Miller (1976) calculate a value of 1. 9-2.0 x 10 ^^gs’^ be- 
tween 8-9 km for a northeast Colorado storm calculated using vertical velocity 
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inferred from dual-Doppler radar data. Auer and Marwitz (1968) present results of 
the cloud base mass flux into 18 hailstorms on the high plains deduced from air- 
craft measurements. Their average value is 2.3 x lO^^gs"^. Therefore, it appears 
that the upward flow deduced from the satellite observations is of a reasonable 
magnitude when compared to calculations and observations on approximately the 
same scale. Inferences about the magnitude of the maximum updraft are presented 
in Section 5. 

4. EXAMPLE OF INTENSE THUNDERSTORM 

On April 24, 1975 a severe thunderstorm complex developed over extreme 
northeastern Oklahoma in the late afternoon and moved into southwestern Mis- 
souri around sunset. The most significant severe weather associated with the sys- 
tem was the Neosho, Mo. tornado which touched down at approximately 0040 
GMT, April 25. By following the evolution of the storm system backwards in time, 
the initial intense convection can be detected and its associated rapid cloud top 
growth calculated. 

The Neosho cloud system was designated cloud 18 as part of a larger study 
of this day. Figure 4 exhibits minimum cloud top Tgg as a function of time for 
cloud 18 in its early stages. The temperature drops precipitously between 2200 
and 2220GMT with a maximum calculated rate of 4K min"E The drop in tem- 
perature between 260 K and 220 K takes only a little more than 15 minutes. This 
type of rapid change emphasizes the importance of short interval data to study 
and monitor thunderstorm activity. 

Using a lapse rate varying from 7.8km'^ at 260K, to 8.6Kkm‘^ at 240K, to 
8.0Kkm‘^ at 21 5 K the vertical velocity was calculated using Equation (1). The 
maximum calculated w is 7.8ms* ^ at about 9 km. 

5. INTERPRETATION OF CALCULATED VERTICAL VELOCITIES IN TERMS 
OF MAXIMUM UPDRAFT 

The vertical velocities presented in the previous sections of this paper are mean 
velocities over an area equivalent to the satellite instantaneous field of view (IFOV), 
which in this case is approximately lOOkm^. In the temperature range 23 5-240 K 
(~8.7km) the 19 observed w’s based on the satellite data ranged from 1.2 to 7.8 
ms"^ with a mean of 3.5ms’C Although these are very large values when com- 
pared to typical synoptic-scale w’s, they are small when compared to maximum 
thunderstorm updraft magnitudes. Thunderstorm updrafts can reach magnitudes 
of 10ms“* very easily and are typically 30 ms'^ in supercell thunderstorms (Brown- 
ing, 1977; Davies-Jones, 1974). These large updraft values are probably representa- 
tive of an area approximately 1 km^. Thus there are two orders of magnitude dif- 
ference in the area covered by the estimated w’s obtained from the satellite data 
in this study and the area covered by the peak updraft velocity. 

Assuming axial symmetry and a knowledge of the shape and size of the radial 
profile of vertical velocity, one can make an estimate of the maximum updraft mag- 
nitude. The formula chosen in this study, based on Kyle et al., (1976) and Schlicht- 
ing (1978) is 


w = Wq (3) 

where w^ is the peak w, r is the radial distance, R is the radius of the updraft, and 
the constant a = 2.3. 
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Integrating Equation (3) over a circular area of radius rj, and dividing the re- 
suit by the area of the integration produces an expression for the mean w over the 
area, i.e., 


- Wo /R\2 

w= — ^—\ [1 - e-^(ri/R)2] ^4^ 

The size of thunderstorm updrafts is highly variable (Browning, 1977). Based 
on the discussion of updraft sizes by Browning (1977) and an examination of cross- 
sections by Kropfli and Miller (1976) and Ray (1976) of thunderstorm vertical mo- 
tions deduced from Doppler radar data, an updraft radius of 3km is reasonable. 

For an area of lOOkm^ (approximately equal to the IFOV), the equivalent 
radius of integration is 5.6km. Thus with R = 3km and rj = 5.6km, 

w=0.12wq. (5) 

Thus, with all the assumptions as to profile shape and updraft size, Equation (5) 
indicates that the mean w of 3.5ms‘i is equivalent to a w^ of 29ms'i and the 7.8 
ms~^ value from the Neosho storm (Section 4) is equivalent to a Wq of 65ms“T 

For an R of 3.5 km, instead of 3km, the 65 ms"^ value in the largest paragraph 
would drop to 46ms"^. It is obvious that the calculated Wq is sensitive to the size 
of the updraft. 

Despite the variability and sensitivity of the w^ calculations, it is evident that 
the satellite observations on a scale of 10km are producing w of up to approxi- 
mately 8ms’^, and this can be interpreted as being roughly equivalent to a maxi- 
mum updraft of 50ms‘i in intense thunderstorms. 

6. SUMMARY 

Rapid-scan (5 minute interval) SMS/GOES IR data have been used to estimate 
thunderstorm top ascent rates for severe and non-severe thunderstorms on two case 
study days. On both days examined (May 6 and April 24, 1975) the thunderstorm 
elements with associated severe weather reports have larger average w’s. The severe 
and non-severe elements had mean w’s of 4.6 and 2.5 ms“^ Intensity of convection 
appears to be correlated with the occurrence of severe weather, and the satellite 
data appear to be capable of quantifying the convection intensity. 

The calculated vertical velocities are representative of an area (100 km^) roughly 
equivalent to the satellite instantaneous field of view (IFOV). Mass flux estimates 
of approximately 2 x lO^^gs^^ are calculated, which are reasonable in comparison 
with other estimates. 

Calculations were performed to estimate the peak updraft velocity from the 
satellite based values (averages over 1 00 km^ areasX With a reasonable value of R 
(3-3.5km), the w of 7.8ms"^ for the Neosho storm produces an estimate of approx- 
imately 50ms’i for Wq. 
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Fig. 1 “Composite vertical velocity 
profiles for May 6, 1975. Numbers in 
parentheses are the number of cases. 



Fig. 3“Frequency distributions of 
estimated vertical velocities at 8.7 km 
for clouds on April 24 and May 6, 
1975. Hatched portion of histogram 
indicates thunderstorms with accom- 
panying severe weather reports. 



Fig. 2-Composite vertical velocity 
profiles for April 24, 1975. 



Fig. 4-Minimum equivalent blackbody 
temperature (Tgg) as a function of 
time for cloud 18 on April 24, 1975. 
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